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Abstract 

Planned originally as a review of CP violation (CPV) in B-mesons 
which covered recent B-factories results these lectures turned out to 
be a bit wider. It is not natural to be limited by CPV in decays and 
mixings of B-mesons and not to speak about the analogous phenomena 
in K-mesons since it is very useful and interesting to study what is 
common and what is different in these systems and why. CKM matrix 
elements are extracted from K and B mixings and decays and the 
deviation from unitarity of CKM matrix may become the place in 
which New Physics will show up. So we discuss this simple and elegant 
piece of Standard Model as well. 

In order to follow these lectures you should be able to write the 
Lagrangian and to draw the Feynman diagrams in a Standard Model 
and to calculate the corresponding amplitudes. 
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1 "The road map 



In Fig.l you can see a set of bounds on the parameters p and fj (p and fj are 
defined by eq.(3.10)) of the Cabibbo-Kobayashi-Maskawa [1, 2] quark mixing 
matrix (CKM). They comprise three circles, two branches of a hyperbola, and 
two straight lines. Three circles originate from the V u b (for definition of V u b 
see eq.(3.4)) measurement (the green one with the center at p = fj = 0, 
see eq.(4.7)), the measurement of Am Bd - the mass difference of two energy 
eigenstates of B d and B& mesons (the red one, eq. (11.16) with the center at 
p — 1, fj — 0) and from the lower bound on Am^ which is the same for 
B s and B s mesons (the yellow one, eq.(11.20)). The hyperbola originates 
from the measurement of CP violation - CPV - in the mixing of i^-mesons 
(see eq.(9.9)). Straight lines come from the measurement of CP asymmetry 
in B d (B° d ) -> J/^K decays (see eq.(14.11)). The fact that all three circles, 
hyperbola and straight lines intersect at one and the same place means the 
triumph of Standard Model. Our aim is to explain in these lectures where 
all these bounds come from. 



2 CKM matrix — where from? 

In constructing the Standard Model Lagrangian the basic ingredients are 
1. gauge group, 2. particle content and 3. renormalizability of the theory. 
There is no such a building block in Standard Model as CKM matrix in 
charged currents quark interactions. CKM matrix originates from Higgs 
field interactions with quarks. The piece of the Lagrangian from which the 
up quarks get their masses looks like: 

A£ up = f^Q l ' L u k R H + c.c. , i, k = 1, 2, 3 , (2.1) 



where 



Ql 






Ur = u' r , u 2 r = c' r , u 3 R =t' R (2.2) 
and H is the Higgs doublet: 



H = I") (2.3) 
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Figure 1: The domains at (p,fj) plane allowed at la from V u t, Am Bd , Ek 
and sin 2f3 measurements. 95%C.L. upper bound from the search of Am# s is 
shown as well. 



The piece of the Lagrangian which is responsible for the down quark 
masses looks the same way: 

AC down = f^Qid k R H + c.c. , (2.4) 

where 

= d' R , d% = s' R , d 3 R = b' R and H a = e ab H^ , (2.5) 
-( M 

£ab ~ { -l o ) ■ 

After SU(2) x U(l) symmetry breaking by the Higgs field expectation 
value < H° >= v from formulas (2.1) and (2.4) two mass matrices emerge: 

M%uiu k R + Mil wn did k R + c.c. (2.6) 



The matrices M up and Md OW n are arbitrary 3x3 matrices; their matrix 
elements are complex numbers. According to the well-known theorem an 
arbitrary matrix can be written as a product of the hermitian and unitary 
matrices: 

M = UH , where H = H + , and UU + = 1 , (2.7) 

(do not mix the hermitian matrix H with the Higgs field doublet) which is 
analogous to the following representation of an arbitrary complex number: 

a = e i4, \a\ . (2.8) 

From eq. (2.7) it is evident that matrix M can be diagonalized by 2 different 
unitary matrices acting from left and right: 

(m u \ 

m c , (2.9) 

m t J 

where rrii are real numbers (if matrix M is hermitian (M = M + ) then we will 
get Ul — Ur). Having these formulas in mind let us rewrite the up-quarks 
mass term from eq. (2.6): 

u % ' L M ik u k R + c.c. = u' L U£U L MU R ~U R u' R + c.c. = u L M diag u R ~ + c.c. = uM diag u , 

(2.10) 

where we introduce the fields Ul and Ur according to the following formulas: 

u L = U L u' L , ur = U r u' r . (2.11) 

Applying the same procedure to matrix Afd own we observe that it becomes 
diagonal as well in the rotated basis: 

d L = D L d' L , d R = D R d' R . (2.12) 

Thus we start from the primed quark fields and get that they should 
be rotated by 4 unitary matrices Ul, Ur, D l and D R in order to obtain 
unprimed fields with diagonal masses. Since kinetic energies and interactions 
with the vector fields A 3 ^, and gluons are diagonal in the quark fields, these 
terms remain diagonal in a new unprimed basis. The only term in the SM 
Lagrangian where matrices U and D show up is charged current interaction 
with the emission of W-boson: 

AC = gW+u' Ll ,d' L = gW;u Ll ^UtD L d L , (2.13) 

and the unitary matrix V = U^D L is called Cabibbo-Kobayashi-Maskawa 
quark mixing matrix. 



5 



3 CKM matrix: angles, phases, parameteriza- 
tion, unitarity triangles 

One can easily check that n x n unitary matrix has n 2 /2 complex or n 2 real 
parameters. The orthogonal n x n matrix is specified by n(n — l)/2 angles 
(3 Euler angles in case of 0(3)). That is why the parameters of the unitary 
matrix are divided between phases and angles according to the following 
relation: 

n 2 _ n(ra-l) _|_ n(n+l) 

(3.1) 

angles phases 

The next question arises: are all these phases physical observables or, in 
other words, can they be measured experimentally. And the answer is "no" 
since we can perform phase rotations of quark fields (u — > e lC -u, d — > e^d 
...) removing in this way 2n — 1 phases of the CKM matrix. The number 
of unphysical phases equals that of up and down quark fields minus one 
since the simultaneous rotation of all up-quarks on one and the same phase 
multiplies by (minus) this phase all the matrix elements of matrix V. The 
rotation of all down-quark fields on one and the same phase acts on V in the 
same way. That is why the number of the "unremovable" phases of matrix 
V is diminished by the number of possible rotations of up and down quarks 
minus one. 

Finally for the number of observable phases we get: 

"("+ 1 )_ (2 „_ 1)= ("- 1 ) 2 ("- 2 ) . p.2) 

As you see for the first time one observable phase arrives in the case of 3 
quark-lepton generations. 

Now a bit of history. At the time when Cabibbo used the mixing angle 
6 C [1] only three quarks (u, d and s) were known, and his suggestion was to 
mix d- and s-quarks in the expression for the charged quark current: 

J£ = ^^(1 +7 5 )[dcos0 c + ssin0 c ] . (3.3) 

In this way he related the suppression of the strange particles weak inter- 
actions to the smallness of angle 9 C , sin 2 6 C pa 0.05. After the establishment of 
Glashow-Salam- Weinberg theory of the weak interactions through the weak 
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neutral current discovery in 1973 and discovery of a charm quark in 1974 it 
became clear, that 2 quark- lepton generations exist. From the point of view 
of GSW theory the mixing of quark generations should be described by the 
unitary 2x2 matrix which according to eqs. (3.1, 3.2) has one angle and zero 
observable phases. This angle is Cabbibo angle. However even before the 
c-quark discovery Kobayashi and Maskawa noticed that in order to describe 
CP-violation (CPV) Standard Model needs at least 3 quark-lepton genera- 
tions since for the first time the observable phase shows up for n = 3 [2]. At 
that time CPV was known only in i^°-mesons and to test KM mechanism 
one needed other systems. Finally almost 30 years after KM model of CP 
violation was suggested it was confirmed by the magnitude of CPV in neutral 
.B-mesons. 

In the case of three generations the matrix of charged currents looks like: 



(uct)i 



v ud 


v us 


v ub \ 




v cd 


v cs 


v cb 


(; 


V td 


V ts 


v tb J 





(3.4) 



where the matrix elements V^ depend on four parameters: three angles and 
one phase. Let us present one possible way of the parametrization of matrix 
V, which is called "standard" parametrization (straightforwardly generaliz- 
able for n > 3 [3]). It is achieved by consequent rotations in planes (12), 
(13) and (23) and the rotation in plane (13) is accomplished by the phase 
rotation. Performed in such a way the phase rotation cannot be removed by 
U(l) rotations of the quark fields: 

(3.5) 



Rl3 — 

and, finally: 
/ 

V = - 

V 



Cl3 





C13C12 

-C23S12 - S23S 13 C 12 e tS 
iS 

s 12 s 23 — c 12 c 23 s 13 e 




(3.6) 



C13S12 s 13 e 
C23C12 - s 12 s n s 23 e l5 S23C13 

— S23C12 — C23Si 3 Si 2 e tS C23C13 



iS 



(3.7) 
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As the next step we can take the experimental values for Vik and extract 
three angles 6>i 2 , O23 , #13 and phase 5 from them. The system is overcon- 
strained; we have more than 4 experimental numbers (see below) and in this 
way one can check how good CKM model in describing data is. However 
it appeared to be useful to reparametrize with the help of the so-called 
Wolfenstein parametrization. Let us introduce new parameters A, A, p and 
i] according to the following definitions: 

A = S12 , A=— , p — coso , 

S12 S12S23 

i] = — ^— sin <5 , (3.8) 

S12S23 

and get expressions for through A, A, p and rj: 



V = 



where 



v ud 


V us 


V u b 


v cd 


v cs 


v cb 


v td 







( l-A 2 /2 A A\ 3 {p-tri)\ 

-\-iA 2 \ 5 f] l-A 2 /2 A\ 2 
V AA 3 (1- p-ifj) —A\ 2 — iA\ 4 fj 1 J 

(3.9) 

Obtaining the last expression we take into account the following hierarchy 
of angles %: s 13 <^ S23 *C s± 2 1 (see below). This last form of CKM 
matrix is very convenient for qualitative estimates and numerically is rather 
accurate. 

The unitarity of the matrix V leads to the following six equations that 
can be drawn as triangles on a complex plane (under each term in these 
equations the power of A entering it, is shown): 

v: d v us + v; d v cs + v&Vu =0 , (3 n) 



~ A ~ A ~ A 1 " 

KdVub + v: d v cb + v t * d v tb = , 

~ A 3 ~ A 3 ~ A 3 

v: s v ub + v*v cb + v*v tb = , 

~ A 4 ~ A 2 ~ A 2 



(3.12) 
(3.13) 
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(3.14) 



V ud V t * d + V US V* + v ub v t t 
~ A 3 ~ A 3 ~ A 3 







(3.15) 



V cd V t * d + V CS V* + V cb V t * b 
~ A 4 ~ A 2 ~ A 2 







(3.16) 



Among these triangles the four are almost degenerate: one side is much 
shorter than two others, and two triangles, expressed by equations (3.12) and 
(3.15), have all three sides of more or less equal lengths, of the order of A 3 . 
These two nondegenerate triangles almost coincide. To prove this statement 
let us note that with very good accuracy V ud = V tb = 1, that is why two 
sides of these triangles have equal lengths and directions. That is why the 
triangles coincide and their third sides should also be equal: 



Now, since V us and —V cd are almost equal to each other (and to A) we 
come to the following result: 



the validity of which can be checked by eq.(3.9). 

So, as a result we have only one nondegenerate unitarity triangle; it is 
usually described by a complex conjugate of our equation (3.12): 



and it is shown in Fig. 2. It has the angles which are called (3, a and 7 
(according to BaBar collaboration) or 1; 2 and 3 (according to Belle 
collaboration). They are determined from CPV asymmetries in 5-mesons 
decays. 

Looking at Figure 2 one can easily obtain the following formulas: 



v; d v cb = v us v* . 



(3.17) 



(3.18) 



v ud v: b + V cd V* b + V td V t 







(3.19) 



(3 = 7r — arg 



a = arg 




(3.20) 



(3.21) 
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V V 

cd cb 

Figure 2: Unitarity triangle 



7 = arg^-^ = 3 . (3.22) 

~ V cb V cd 

Angle (3 (0i) was directly measured through time dependent CPV asym- 
metry in B,i — > J/^/K decays, a (02) has been measured recently with a 
rather poor accuracy from CPV asymmetries in B d — > -n + it~ decays and, 
finally, B s decays could be important to determine angle 7 (0 3 ). 

Let us make three final remarks about the unitarity triangle: 

1. In standard parametrization which we use V c dV* b is almost real; 

2. Multiplication of any up quark field on a phase does not change the uni- 
tarity triangle while multiplying d- or 6-quark field on a phase we rotate 
it as a whole not changing its angles which are physical observables; 

3. Usually a rescaled triangle is used. We get it dividing all three sides 
by IV^V^I « A\ 3 . In this way the length of the triangle basement 
becomes equal to one while two other sides have the length of the order 
of one. 

Four quantities are needed to specify CKM matrix: s i2 , S13, s 23 and 5, or 
A, A, p, i]. Knowing more we are checking Standard Model and looking for 
New Physics. 
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4 V us , V ch , V uh - first circle 

The most precise value for the quantity V us follows from the extrapolation of 
the formfactor of K — > nev decay f+(q 2 ) to the point q 2 = 0, where q is the 
lepton pair momentum. Due to the Ademollo-Gatto theorem corrections to 
the CVC value /+(0) = 1 are of the second order of flavor SU(3) violation, 
and these small terms were calculated in [4]. As a result of this analysis 
Particle Data Group gives the following value [5] : 

V us = A = 0.2196 ± 0.0026 . (4.1) 

The accuracy of A is high: the other parameters of CKM matrix are 
known much worse. 

The value of V cb is determined from the inclusive and exclusive semilep- 
tonic decays of 5-mesons with charmed particles production. At the level of 
quarks b — > ceu transition is responsible for these decays. 

According to PDG [5]: 

V cb = (41.2±2.0)10^ 3 ; (4.2) 

and the error is dominated by a theoretical one. 

From eq.(3.9) with the help of (4.2) for parameter A we get: 

A = = 0.85 ±0.04 . (4.3) 
A 

From the formula for the semileptonic width: 

Tsl~ \V cb \ 2 {m b -m c f (4.4) 

it is clearly seen that 4% error in V cb corresponds to the knowledge of quark 
mass difference with the very high accuracy: 

m h -m c = 3 GeV±60 MeV . (4.5) 

Vast literature is devoted to the different aspects of the V cb determination. 
Not going into details let us only note that at the moment the error in V cb 
is not a "bottle neck" since the errors in other relevant quantities (V ub , Vu, 
sin 2/3) are significantly larger. 



11 



The value of \V u b\ is extracted from the semileptonic 5-mesons decays 
without the charmed particles in the final state which originated from b — > 
ulv transition, b — > civ decays are approximately 100 times more probable 
than b — > ulv and for their suppression the high energy charged lepton tail is 
examined (the energetic leptons cannot accompany heavy D-mesons). The- 
oretical analysis of such semiinclusive decays is highly involved and it leads 
to a large theoretical uncertainty [5]: 

| ^j- |= 0.40 ±0.08 . (4.6) 

There are hopes to considerably diminish this error using the exclusive 
modes B — > nev and B — > uev. (The recent study of the B — > ixcv and 
B — > pev decays by CLEO collaboration gives 10% smaller central value of 
V u b with practically the same error [6]. 

Using equation (3.9) we obtain a bound on the parameters p and fj: 

^p2 + f} 2 = o.4o ± 0.08 , (4.7) 

which produces a circle on (p, fj) plane with the center at the point (0,0). 
The area between such two circles in Fig.l corresponds to the domain allowed 
at one sigma. 

The detailed articles on V u j, and V c b determination with bibliography can 
be found in [5]; for the recent review see [7]. 

To finish this section let us look at the value of | V u d | . Its value extracted 
from the neutron decays has the smallest theoretical uncertainty [8]: 

| V ud |= 0.9713(13) ; (4.8) 

using this value as well as (4.1), (4.2) and (4.6) we obtain: 

I Vud | 2 + | V us | 2 + | V ub | 2 = 0.9434(25) + 0.0482(11) + 0.00001 = 

= 0.9916(27) , (4.9) 

and the unitarity of CKM matrix is violated at 3a level. This can be a signal 
for New Physics, while it may be a statistical fluctuation as well. PDG 
discussing the value of V u d takes into account the nuclear beta decays as well 
resulting in [5]: 

| V ud |= 0.9734(8) , (4.10) 
which diminishes the violation of unitarity to 2a. 
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5 CPV: history; why phases are relevant 



In 1956 Lee and Yang in order to solve 6 — r problem suggested that P-parity 
is broken in weak interactions [9] . Soon it was noted that C-invariance is also 
broken in a proposed theory [10]. 

In 1957 looking for a way to resurrect P-invariance L.D. Landau stated 
that weak interactions should be invariant under the product of P reflection 
and C conjugation. He called this product the combined inversion and ac- 
cording to him it should substitute P-inversion broken in weak interactions 
[11]. In this way the theory should be invariant when together with changing 
sign of the coordinate, f — > — f, one changes an electron to positron, proton 
to antiproton and so on. 

It is clearly seen from paper [11] that according to Landau CP-invariance 
should become a basic symmetry for physics in general and weak interactions 
in particular. 

This beautiful picture did not stop experimentalists (or may be even 
stimulated) and in 1964 CP violating decay of the long-lived neutral kaon on 
two pions was discovered [12]. The question of the violation of P-parity in 
weak interactions which stimulated Landau was: why is P violated? Landau's 
answer to the question "Why is parity violated in weak interactions" was: 
because CP, not P is the fundamental symmetry of nature. A modern answer 
to the same question is: because in P-invariant theory with the Dirac fermions 
the gauge invariant mass terms can be written for quarks and leptons which 
are not protected of being of the order of Mqut or Mpi anc k. So in order to 
have our world made from light particles P-parity should be violated. 

Kl — > 2ix decay discovered in 1964 occurs due to CPV in the mixing of 
neutral kaons (e ^ 0). Only thirty years later the second major step was 
done: direct CPV was observed in kaon decays [13]: 

T(K L -» vr+vr-) T(K L -> ttV) 

T(K S -> 7T+7T-) T T(K S -> 7T°7T°) ' T K > 

Finally, in the year 2001 CPV was for the first time observed out of the 
decays of neutral kaons: the time dependent CP-violating asymmetry in B° 
decays was measured [14]: 

n(f) _ N(B° - J/*K S(L) ) - N(B° -» J/*K S(L) ) 
^ > ^(fio J/mK s(L) ) + N{B* - J/*K S{L) ) T ' 1 • ] 
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Starting from the year 1964 we know that there is no symmetry between 
particles and antiparticles. In particular, the C-conjugated partial widths 
are different: 

F(A -> BC) ^ T(A -> BC) . (5.3) 

However CPT (deduced from the invariance of the theory under 4-dimen- 
sional rotations) remains intact. That is why the total widths as well as the 
masses of particles and antiparticles are equal: 

M A = M A , T A = T A (CPT) . (5.4) 

The consequences of CPV can be divided into macroscopic and micro- 
scopic. CPV is one of the three famous Sakharov's necessary conditions to 
get as a result of evolution of charge symmetric Universe a charge nonsym- 
metric one [15]. In these lectures we will not discuss this very interesting and 
well developed branch of physics, but will deal with CPV in particle physics 
where the data obtained up to now confirm Kobayashi-Maskawa model of 
CPV. New data which should become available in coming years could well 
disprove it clearly demonstrating the necessity of physics beyond the Stan- 
dard Model. Let us note that there exist many excellent reviews on CPV; a 
partial list of them can be found in ref. [16], which starts from the lectures 
on this topic given at ITEP Winter Schools in a chronological order. 

The next question we wish to discuss is why the phases are relevant for 
CPV. Let us take the simplest example of one scalar particle decaying in two 
scalar particles. The interaction Lagrangian is: 

£ = \AB*C* + \*A*BC . (5.5) 

Performing CP conjugation (changing particles to antiparticles and re- 
flecting space coordinates) we obtain: 

£ C p = XA*BC + \*AB*C* , (5.6) 

which coincides with the original Lagrangian only if A is real, A* = A. That 
is why the complex coupling constants are necessary for CPV. However, this 
complexity is not enough: in our example we can redefine field A by phase 
rotation, A = e l ^A', in this way making A real. However, adding one new 
field (say D) we introduce many new couplings (\iDBC + \2DAB + \ 3 DAC) 
while by rotation D = e ta D' a phase of only one coupling constant can be 
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eliminated. Thus to get CPV several fields are needed. The last statement 
is illustrated in Standard Model by the fact that at least three generations 
are needed to get CPV through the phases in the quark mixing matrix. 

6 M° — M° mixing; CPV in mixing 

In this part the general formulas for the meson-antimeson mixing will be 
derived. In order to mix the mesons must be neutral, not coincide with its 
antiparticle and decay due to the weak interactions. There are four such 
pairs: 

K°(sd) - K°(sd) , D\cu) - D\cu) , 

B° d (bd) - B° d (bd) and B° s (bs) - B° s (bs) . 

Fast t — > bW decay prevents forming t-quark containing hadrons. 
Mixing occurs in the second order in weak interactions through the box 
diagram which is shown in Fig. 3 for K° — K° pair. 



S U,C,t d 




Figure 3: K — K mixing 



The effective 2x2 hamiltonian H is used to describe the meson-antimeson 
mixing. It is most easily written in the following basis: M° = ( ?. ] , M° = 




. The meson-antimeson system evolves according to the Shroedinger 



equation with this effective hamiltonian which is not hermitian since it takes 
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meson decays into account. So, H = M — |T, where both M and T are 
hermitian. 

According to CPT invariance the diagonal elements of H are equal: 



< M° | H | M° >=< M° | H | M° > 
Substituting into the Shroedinger equation 



tp - function in the following form: 



1> 



we come to the following equation: 



q 6 ' lXt 



I 


M — 


ip 

2 L 


M 12 - 


i p 


\ 


/ P \ 




( p \ 
















= X 




\ 


M* 2 - 


i p* 
2 1 12 


M — 


ir 

2 


/ 


V 9 J 




V i J 



from which for eigenvalues (A±) and eigenvectors (M ± ) we obtain: 



M* _ IF* 
12 2 12 

'12 - ^ I 12 



(6.1) 



(6.2) 



(6.3) 



(6.4) 



A± = M - l -Y ± J (M 12 - l -Y l2 ){Mt 2 - % -Y\ 2 ) , (6.5) 



(6.6) 



\ M + = pM° + qM° q _ 
\ M_= pM° - qM° ' p ~ \\ Mi; 

Multiplying M° by e i<?i we are changing the phase of the ratio q/p, that 
is why arg ( - ) is not a physical observable. 



CP transforms field M° in the following way: 

CPM° = e ia M° , 



(6.7) 



being defined up to the arbitrary phase a. We can rotate field M° by phase 
a/2, removing phase a from the definition of CP transformation. In this way 
we come to the standard definition of CP transformation: 



CPM° = M° , 



(6.8) 
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simultaneously loosing freedom of the arbitrary phase rotation of field M°. 
The eigenstates M + and M_ are not orthogonal in general case: 



< M+ | M_ >=| p | 2 - 
However if there is no CPV in mixing, then: 



?lVo 



< M° I if I M u >=< M u I ff I M u > 



rO 



o 



rO 



m 12 - -r 12 



JKi 12 2 12 



that is why: 



p 



1 , < Mi I M_ >= 



(6.9) 

(6.10) 
(6.11) 

(6.12) 



We observe the one-to-one correspondence between CPV in mixing and 
nonorthogonality of the eigenstates M + and M_. According to Quantum 
Mechanics if two hermitian matrices M and V commute, they have a common 
orthonormal basis. Let us calculate the commutator of M and T: 



M 12 r* 2 - MUT 



[M,Y\ = 



12 1 12 











M* 2 r 12 - M 12 r* 12 J 



(6.13) 



It equals zero if the phases of Mi 2 and Ti 2 coincide modulo it. So, for 
[Mr] = we get | q/p \= 1, < M+ | M_ >= and there is no CPV in the 
meson-antimeson mixing. And vice versa. 

Introducing quantity e according to the following definition: 

q 1 — e 

p 1 + e 

we see that if Re e ^ 0, then CP is violated. From (6.6) for the eigenstates 
we obtain: 



(6.14) 



Ma. 



M_ 



1 



v/iT 



e I 2 



M° + M° M° - M° 
7= +^ 

M°-M° 



V2 



+ e- 



V2 
M° + M° 



(6.15) 



If CP is conserved, then Re e = 0, M + is CP even and M_ is CP odd 
(CP transformation should be defined according to (6.7) - (6.8)). If CP is 
violated in mixing, then Re e ^ and M + and M_ get admixtures of the 
opposite CP parity and become nonorthogonal. 
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7 Space-time pattern of K° — K° 
(B° — B° , — v e ) oscillations 

Neglecting CPV for the states with the definite masses and width we obtain: 

= £(*° + *°) • (71) 

where CP|^° >= \K° >, CP\K° >= \K° >. 

A plain wave which describes the propagation of the superposition of Ki 
and K 2 with definite energies and momenta is: 

V(x,t) = e~ iElt+ipiX | K x > +e - lE * t+ip2X | K 2 > , (7.2) 

7-i2 2 2 77i2 2 2 

E 1 -p 1 = m 1 , E 2 -p 2 = m 2 , 

and we consider the propagation in the direction of axis x. We should impose 
a boundary condition stating that at point x = continuously in time K° 
is produced (or B°, or v e ). The only way to get this is by putting E 1 = E 2 , 
having at x = 0: 

e -iEt | Ri > +e -iEt \ K2>= e -iB*[| > + | X 2 >] = e - iE * | > . 

(7-3) 

If we do not impose this boundary condition and take E 1 ^ E 2 , then after 
time r ~ e 1 -e 2 a ^ ^ e P ^ 11 ^ ^ = -ft' will emerge, while in the experiment 
only the events when K° at x = is produced are selected. For example, in 
CPLEAR we have: pp — > nK~K° and K° is tagged by K~ , showing that 
K° was produced in the interaction point, and not K°. Another example: in 
the reaction pn — > pAK° together with A hyperon K° is produced and not 
K° because of strangeness conservation in strong interactions. In the case of 
antineutrino production in nuclear reactor initiated by n — > pe~v e transition 
it is always electron antineutrino and never muon antineutrino. 

Substituting E 1 = E 2 = E in equation (7.2) we get: 

V(x,t) = e - iEt+ipiX [\ K 1 > + e iip2 - pi)x | K 2 >] . (7.4) 
For the phase factor in brackets we have: 



Vi~V\ — v E 2 — m 2 — J E 2 — ml 
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= tJe 2 -ml + ml-ml- ^Je 2 -m\ = m \ ™" 1 , (7.5) 

2p 

where we take into account that for mesons (as well as for neutrinos) always 
E 2 — m\^> m\ — m\. Therefore for a probability to detect K° at a distance 
x from a production point we have: 



P K o K o — - 



1 + cos( 



m\ 



m. 



2p 



x) 



Am<gCr; 



Am 
1 + cosi— — x] 



Pi 



(7.6) 



j3 = v/c, 7 = 1/^1-^ , 



where the second equality in eq.(7.6) holds for mesons, but (may be) not for 
neutrinos. 

There is no surprise that both neutrino and meson oscillations are de- 
scribed by the identical formulas. In the first paper where the neutrino 
oscillations were considered [17], B.M. Pontecorvo did it in full analogy with 
K° — K° oscillation analysis of Gell-Mann and Pais. Since at that time the 
second (muon) neutrino was not yet discovered Pontecorvo considered v e — v e 
oscillations which were allowed because V — A theory was not established yet 
(in V — A theory v — v oscillations are forbidden by chirality conservation). 
In this particular case the diagonal elements of v — v mixing matrix should be 
equal due to CPT and mixing is maximal, Q v = 7r/4, just as for K°(B d ,B s ) 
mesons. 



8 K° - K° mixing, Am LS 

Ti2 for the K° — K° system is given by the diagram shown in Fig. 4. With 
our choice of CKM matrix V us and V U( i are real, so Ti2 is real. 

Myi is given by dispersion part of the diagram shown in Fig.5. Now all 
three up quarks should be taken into account. 

To calculate this diagram it is convenient to implement Glashow-Illiopulos- 
Maiani compensation mechanism from the very beginning, subtracting zero 
from the sum of the fermion propagators: 

V V- V* 

V V* V V* \L V* Z- 1 ls id 

V us V u d _|_ V cs V c d _|_ Vjs Vjd _ J. -|\ 

p - m u p - m c p - m t p 
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u 

* Xr 



W ( ( W 



-K — >- 



d us 

Figure 4: Quark diagram responsible for T 12 in K° — K° system 

s u,c,t d 
* 



W ( ( W 



d u,c,t s 

Figure 5: Quark diagram responsible for M 12 in K° — K° system 

Since w-quark is with good accuracy massless, m u rs 0, its propagator 
drops out and we are left with the modified c- and t-quark propagators: 



P ~ m c>t {p 2 - m 2 t )p 



(8.2) 



The modified fermion propagators decrease in ultraviolet so rapidly that 
one can calculate the box diagrams in the unitary gauge, where W-boson 
propagator does not decrease (one can demonstrate that the diagrams with 
the charged higgs exchange which occur in the renormalizable R gauges, 
and for which eq.(8.1) does not work since the vertex of higgs emission is 
proportional to the quark masses, becomes zero in the limit £ — > oo, which 
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corresponds to the unitary gauge). 

We easily get the following estimates for three remaining diagram contri- 
butions in M 12 : 

(cc) : X 2 (l + 2ifjA 2 X 4 )G 2 F m 2 c , (8.3) 



(ct): X^l-p-t^GlmlH^) 2 , (8.4) 

m c 

(tt) : X w (l-p-tf]) 2 G 2 F m 2 . (8.5) 

Since m c w 1.3 GeV and m t 175 GeV we observe that the cc diagram 

dominates in i?eM 12 while imM 12 is dominated by (tt) diagram. The real 
part dominates in M 12 : 

A 8 f^V^O.l . (8.6) 



ReM 12 \m c J 
The explicit calculation of the cc exchange diagram gives: 

££=2 = -^r^;(^(i + 7.)d) WO£ , (8.7) 

where g is SU(2) gauge coupling constant, g 2 /8M^ v = Gf/ V%, and the factor 
771 takes into account the hard gluon exchanges. Since 

M 12 - l -T 12 =< K° I H eff I K° > /2m K l , (8.8) 

(here H e ^ = —C e /£ =2 ) we should calculate the matrix element of the product 
of two V — A quark currents between K° and K° states. Using the vacuum 
insertion we obtain: 

< K° | s 7a (l +l5 )ds la (l + l5 )d | 7?° >= 



1 Factor l/2mjf appears when taking a square root from the quadratic in meson masses 
mesonic Hamiltonian in order to obtain linear in meson masses Hamiltonian which enters 
the Shroedinger equation: 

(m - §r) 2 < if 1 H e ff \k°> ] 1/2 _ f m - |r <g0| ^" |J?0> N 



< | # e ^' | if > (m - ^r) 2 



2m jc 2 



here = -£^/ =2 . 
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= -B K < K° | s 7a (l + ls )d | > x (8.9) 

< | s 7a (l + 75 )^ I ^° >= - , 

where i?x = 1 if the vacuum insertion saturates this matrix element. 
With the help of eq.(6.5) we obtain: 

m s - m L - l -(F s - T L ) = 2[ReM 12 - % -Y l2 ] , (8.10) 

where S and L are the abbreviations for K s and K L , short and long-lived 
neutral i-T-mesons respectively. For the difference of masses from (8.7 - 8.10) 
we get: 

— a G 2 F B K f K m K 2 T/ 2 T/ * 2 / 011 \ 
m L -m s = Am LS = — 2 Vim c V cs V cd . (8.11) 

Constant fx is known from K — > Iv decays, fx = 160 MeV. Gluon 
dressing of the box diagrams in 4 quark model in the leading logarithmic 
(LO) approximation was calculated in [18], r]f° = 0.6. It appears that the 
subleading logarithms are numerically very important [19], ??f LO = 1.3 ±0.2, 
the number which we will use in our estimates. We take = 1 ± 0.1 
assuming that the vacuum insertion is good numerically, though the smaller 
values of B K can be found in literature as well (see [20]). 

Experimentally the difference of masses is: 

Aml x | = 0.5303(9) • 10 10 sec" 1 ; (8.12) 
Substituting the numbers in eq. (8.11) we get: 

A theor 

A^ = °' 6±0 ' 2 ' (8 ' 13) 

and we almost get an experimental number from the short-distance contribu- 
tion described by the box diagram with c-quarks. As V cs and V c d are already 
known nothing new for CKM matrix elements can be extracted from Atjils- 
Concerning the neutral kaon decays we have: 

\ 2 w T 5 = 1.1 • 10 10 sec" 1 (Am L5 ^r s /2) , (8.14) 
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since Tl <C r#, = 2 • 10 7 sec" 1 . Kl is so long-lived because it can decay 
only into 3 particles finite states (neglecting CPV) and for the decays into 3 
pions the energy release is small: 

m K - m 3n = 490 MeV - 3 • f 40 MeV = 70 MeV . (8.15) 

Ks rapidly decays to two pions which have CP= +1. 



9 CPV in K° - K () : K L 2tt , e K -hyperbola 

CPV in K° — K° mixing is proportional to the deviation of | q/p | from one; 
so let us calculate this ratio according to eq. (6.6) taking into account that 
Ti2 is real, while M 12 is mostly real: 

q _ _ ilmMu _ 2iImM 12 _ , . 

p - m 12 - |r 12 ~ m L - m s + |r 5 • 1 • j 

In this way for the quantity e introduced in eq. (6.14) we obtain: 

HmMi2 

Branching of CP-violating K L — > 2n decay equals: 



Br(K L 



2A + Br(Kl _ = £<£l^) = E^EjM = 

J- K L J- K s ^2w J- yJ^L) 

Voo | 2 T(K S -> 27r°)+ | | 2 T(K 5 -> tt+tt-) r(gg) _ 

r(K 5 -> 2tt°) + r(# s -> tt+tt-) r(K L ) ~ 



v/2 



e i?eA 2 V2 



5.17(4) • 10" 8 sec 
0.894(1) • 10- 10 sec 



« 578(1 - 0.08) | e | 2 = 3.02(3) • 10~ 3 

where the last number is the sum of Kl — > 7r + 7r~ and — > 7T 7T branching 
ratios and experimental values of the i^s and widths are used. We also use 
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the approximate relation Sq — S 2 ~ 45° which follows from the analysis of 71 — tt 
scattering. In eq. (9.3) the small factors of the order of e'/e are neglected 
while the enhanced term j^j^j = 22.24 is taken into account (e'/e is almost 
real). This term originates from direct CPV in kaon decays, see the next 
section. We used the following approximate formula: i] 00 = e + i 1 ^^ , and 
estimated the ratio of amplitudes using eq. (10.10), in which only induced by 
QCD penguin term ~ ImA / ReA is taken into account. The electroweak 
penguins partially cancel QCD penguin; omitting them we underestimate 
I £ I a bit. 

In this way the experimental value of | e | is determined, and for our 
theoretical result described by eq. (9.2) we should have: 



I ImM 12 
v^Am 



= 2.38(1) • 10" 3 (9.4) 



LS 

To calculate ImMi 2 one should use eq. (8.8) and the expression for 
£a/=2 which follows from the calculation of the box diagrams. As we already 
demonstrated (tt) box gives the main contribution in /mM 12 . For the first 
time it was calculated explicitly not supposing that m t <C in the paper 
[21] (a year later the same result was independently obtained in [22]): 

ImM 12 = - GlB ^* mK mh^lVtl) x 7(0 , 

r^-ne + 4 j^inn (m t y 

where factor r] 2 which takes into account the gluon exchanges in the box 
diagram with (tt) quarks was found in the same paper [21] in the leading 
logarithmic approximation: 7]^° = 0.6. This factor is not changed substan- 
tially by subleading logs [19]: r) 2 LO = 0.57(1). Concerning factor i] 3 which 
is responsible for gluon dressing of (ct) box and which was as well found in 
[21]: rfe° = 0.4, the subsequent approximation is [19]: LO = 0.47(4). 

Let us present the numerical values for the expression in figure brackets 
in eq. (9.5) for several values of the top quark mass (this factor was obtained 
in paper [21] and later in literature was named Inami-Lim factor, an evident 
example of the famous Arnold principle: "If a notion bears a personal name, 
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then this name is not the name of the discoverer"): 



1, m t = 0, £ = 
{ } = 0.55 , £ = 4.7 , which corresponds to m t 
0.25 , m t = £ = oo 



175 GeV 



(9.6) 



It is clearly seen from (9.5) and (9.6) that the top contribution to the box 
diagram is not decoupled (it does not vanish) in the limit m t — > oo. Three 
years earlier the analogous non-decoupling of t-quark contribution through 

loops to the quantity p = g i j M $ was observed in [23] . 

One can easily get where this enhanced at m t — > oo behaviour originates 
from the estimating box diagram in 't Hooft-Feynman gauge. In the limit 
m t ^> mw the diagram with two charged higgs exchanges shown in Fig. 6 
dominates, since each vertex of higgs boson emission is proportional to m t . 




Figure 6: Box diagram the contribution of which is enhanced as m 2 in the 
limit m t ^> mw 

For the factor which multiplies the four-quark operator from the diagram 
shown in Fig. 6 we get: 



,m 



t\4 



(-) 



p 2 - M w y 



V 



p 2 — m 2 



~ ( ) — n = G F m t ) 

v m 2 t t 1 



(9.7) 



where v is the Higgs boson expectation value. 

Substituting eq. (9.5) into (9.4) and substituting the numbers we obtain: 



Hhcor= 0.007577(1 -p)(l± 0.1) , 



(9.8) 
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fj(l -P) = 0.32(3) , 

where 10% uncertainty in the value of Bk = 1 ± 0.1 dominates in the error. 
Taking into account (ct) and (cc) boxes we get the following equation: 



77(1.52 — p) = 0.32(3) , 



(9.9) 



which determines the hyperbola in Fig. 1. 

Let us make two comments before finishing this section. 

a) Why is Ek small? From eqs. (9.4), (9.5) and (8.11) we obtain the 
following estimate for e K : 



It means that Ek is small not because CKM phase is small, but because 
2x2 part of CKM which describes the mixing of the first two genera- 
tions is almost unitary and the third generation almost decouples. We 
are lucky that the top quark is so heavy; for m t ~ 10 GeV CPV would 
not be discovered up to now (the last statement is true only if the 
mixing angles of the first two generations with the third one are not 
proportional to 1/ ^Jrn~t - in the opposite case Ek would not depend on 
m t and its smallness would not have any qualitative explanation). 

b) If the masses of any 2 up (or down) quarks are equal, then CPV escapes 
from CKM matrix. Indeed, for m c = m t we have: 



ImM 12 ~ Im[V ts V t * d + V cs V c * d ) 2 = Im[V us V: d } 2 = . (9.11) 

10 Direct CPV in K decays, e' / (| | |^ 1) 



Let us consider the neutral kaon decays into two pions. It is convenient to 
deal with the amplitudes of the decays into the states with definite isospin: 



e k ~ 



m 2 X 10 7](l — p) 
m1\ 2 



(9.10) 



A(K° -> 7T+7T-) 




(10.1) 



A(K° -> 7T + 7T-) 




(10.2) 
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A (K° - 7r°7T ) = J^a 2 e~^e^ - -^e^e* , (10.4) 
V 3 v3 

where "2" and "0" are the values of (tttv) isospin, £2,0 are the weak phases 
which originate from CKM matrix and #2,0 are the strong phases of tttt- 
rescattering. If the only quark diagram responsible for K — > 27r decays were 
charged current tree diagram which describes s — > uud transition through 
VT-boson exchange, then the phases would be zero and it would be no CPV 
in the decay amplitudes (the so-called direct CPV). All CPV would originate 
from K° — K° mixing. Such indirect CPV was called superweak. However 
in Standard Model CKM phase penetrates into the amplitudes of K — > 2ir 
decays through the so-called "penguin" diagram shown in Fig. 7, and £0,2 
are nonzero leading to direct CPV as well. 




W 

Figure 7: Penguin diagram responsible for direct CPV in Standard Model 
From eqs. (10.1) and (10.2) we get: 

\/2 

o 

so for direct CPV to occur through the difference of K° and K° widths at 
least two decay amplitudes with different CKM and strong phases should 
exist. 

In the decays of Kl and K$ mesons violation of CP occurs due to that 
in mixing (indirect CPV) and in decay amplitudes of K° and K° (direct 
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CPV). The first effect is taken into account in the expression for Kl and Kg 
eigenvectors through K° and K°. With the help of eq. (6.15) we readily 
obtain: 

K° + K° K°- K° 
Kg = — + e- 



V2 

K°-K° 
K L = = — + e 



y/2 
K° + K° 



(10.5) 



y/2 V2 ' 

where we neglect ~ e 2 terms. For the amplitudes of Kl and Kg decays into 
7r + 7r~ with the help of eqs. (10.1), (10.2) we obtain: 



A{K L 



' +0 = 7=2 



■^=e i52 2isinf 2 + ^=V2e i5o 2ism^ 
v 3 v 3 



+ 



+ 



V2 . 



MKg - 7T+7T-) = -L 



-^=e i& 2cose2 + -^=V2e i<5o 2coseo , (10.6) 
_ v 3 v 3 

^e*2cos6 + ^v^e < *2oos& , (10.7) 

where in the last equation we omit the terms which are proportional to the 
product of two small factors, e and sin£ ,2- For the ratio of these amplitudes 
we get: 

A(K L -> TT+TT") _ .sinfo ie i( - s ' 2 ~ St)) a 2 cos£ 2 [sin^ sin£ 
A(K S — > 7T+7T") cos^o a cos£ 



COS £ 2 COS £ 

(10.8) 



(1 

where we neglect the terms of the order of (a2/ao) 2 sin£o,2 because from 
AT = 1/2 rule in i^-meson decays it is known that 02/00 ~ 1/22. 

The analogous treatment of K L S — > 7r°7r° decay amplitudes leads to: 

sin £ 2 sin £ 



^-je<(fc-*0y/2 a2COs6 
cos^o ao cos Co 



A(K L -> 7r°7r°) _ _ ._ 
7700 ~ A(K S - 7r%°) " £ + 'cose 

The difference of r)± and 7700 is proportional to 



cos { 2 cos £ 



y/2 a cos £ 



sin sin £ 

cos £2 cos £ 



(10.9) 
(10.10) 
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y/2 ReA 



ImA 2 ImAo 



1 __ e i(S2-5 ) 1 



V2 ReA c 



ImA 2 - ^IitlAq 



ReA 2 ReA 

where A 2fi = e^ 2 - a 2 ,o- 

Introducing quantity e according to the standard definition: 

IuiAq 
£ = £ + l R^ ' 

we obtain: 

V+- =e + e' , r l00 =e-2e' . (10.11) 

The double ratio rj±/r] 00 was measured in the experiment and its difference 
from 1 demonstrates direct CPV in kaon decays: 

/ ,\ cxp 

( - J = (1.8±0.4)10~ 3 . (10.12) 

The smallness of this ratio is due to (1) the smallness of the phases pro- 
duced by the penguin diagrams shown in Fig. 7 and (2) smallness of the 
ratio of a 2 /a ~ ReA 2 /ReA . 

Let us estimate the value of e' . The penguin diagram with the gluon 
exchange generates K — > 2n transition with AI = 1/2; those with 7- and 
Z-exchanges contribute to AI = 3/2 transitions as well. The contribution 
of electroweak penguins being smaller by the ratio of squares of coupling 
constants is enhanced by the factor ReA / ReA 2 = 22, see the last part in 
eq. (10.10). As a result a partial compensation of QCD and electroweak 
penguins occur. In order to obtain an order of magnitude estimate let us 
take into account only QCD penguins. From eq. (10.10) we obtain the 
following estimate for the sum of the loops with t- and c-quarks: 

I e' I- J-±k = -i_^Wl n( ^)^A 5 - 
22V2cos£ 22 V2 12?r y m c ' 

^ 10 - 5 ^K) ln( ^ )2 (1013) 

12tt y m c J y ' 

Taking into account that | e |« 2.4 • 10~ 3 we see that the smallness of the 
ratio of e'/e can be readily understood. 

In order to make an accurate calculation of e'/e one should know the 
matrix elements of the quark operators between if-meson and two 7r-mesons 
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which with the present knowledge of low-energy QCD is not possible. That 
is why a horizontal strip in Fig. 1 which should correspond to eq. (10.13) has 
too large width and is not shown. Nevertheless we discussed direct CPV in 
this section since it will be important for 5-mesons. For more details about 
e'/e estimations see review papers [24]. 



11 — and B® — B® mixing — two circles 

The £>-meson semileptonic decays are induced by a semileptonic 6-quark 
decay, b — > l~vc{l~vu). In this way in the decays of B° mesons l~ are 
produced, while in the decays of B° mesons l + are produced. However B° 
and B° are not the mass eigenstates and being produced at t — they start 
to oscillate according to the following formulas: 

B°(t) = - J£ 5° + 11 _1 5° , (11.1) 

B\t) = e - 3° + ^- ^ 5° ; 

for derivation use eq. (6.5 - 6.6). 

That is why in their semileptonic decays the "wrong sign leptons" are 
sometimes produced, l~ in the decays of the particles born as B° and Z + in 
the decays of the particles born as B°. The amount of these "wrong sign" 
events depends on the ratio of the oscillation frequency Am and £>-meson 
lifetime T (unlike the case of K-mesons for 5-mesons Ar <^ T, see below). 
For Am/r » 1 a large number of oscillations occurs, and the number of 
"the wrong sign leptons" equals that of a normal sign. If Am T, then 
.B-mesons decay before they start to oscillate, and the number of "the wrong 
sign events" is power suppressed. The pioneering detection of "the wrong 
sign events" by ARGUS collaboration [25] in 1987 demonstrates that Am 
is of the order of T, which in the framework of Standard Model could be 
understood only if a top quark is unusually heavy, m t > 100 GeV. Fast 
B° — B° oscillations made possible the construction of asymmetric 5-factories 
where CPV in B° decays was observed. (At the end of this introduction let 
us mention that UA1 collaboration saw the events which were interpreted as 
a possible manifestation of B° s — B° s oscillations [26].) 
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Integrating the probabilities of B° decays in l + and I over t, we obtain 
for "the wrong sign lepton" probability: 

N B o^i- x 



W, 



B o_ B o = 



N B o^i- x + N B o^i+ x 



q 12 / Am\2 

7 ' 1 (11.2) 



" 2+ (^p)2 + | 1 |2 (Aj2)2 ' 

where we neglect Ar, the difference of B + - and i?„-mesons lifetimes. Pre- 
cisely according to our discussion for Am/T > 1 we have W = 1/2, while 
for Am/T « lwe have W = l/2(Am/r) 2 (with high accuracy | p/q \— 1, 
see below). 

For B° decays we get the same formula with the interchange of q and p: 

N B o^i+ x 



W, 



lv B o^i+ x + l\ B o^i- x 

I p 12 / Arn \2 

= 1 1 1 ^ r 1 (H 

2 + (^) 2 + | 2 | 2 (^p) 2 ' K ' J 

In ARGUS experiment B- mesons were produced in T(4S) decays: T(4S*) — 
BB. For T resonances J PC = 1 , that is why (pseudo) scalar 5-mesons are 
produced in P-wave. It means that BB wave function is antisymmetric at 
the interchange of B and B. This fact forbids the configurations in which due 
to B — B oscillations both mesons became B, or both became B. However 
after one of the 5-meson decays the flavor of the remaining is tagged, and it 
oscillates according to eqs. (11.1). 

If the first decay is semileptonic with Z + emission indicating that a decay- 
ing particle was B°, then the second particle was initially B°. With the help 
of eqs. (11.2 - 11.3) and taking | p/q \— 1 we get for the relative number of 
the same sign dileptons born in semileptonic decays of £>-mesons, produced 
in T(4S) -> BB decays: 

JW + iW W Am 

N l+l - ~ 1-W ~ 2 + x 2 ' x ~ T ' 1 > 

Let us note that if B° and B° are produced incoherently (say, in hadron 
collisions) a different formula should be used: 

N1+1+ + Ni-i- 2W-2W 2 x 2 (2 + x 2 ) 



Ni+i- 1-2W + 2W 2 2 + 2x 2 + x 4 



(11.5) 
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In the absence of oscillations (x = 0) both (11.4) and (11.5) are zero; for 
high frequency oscillations (x ^> 1) both of them are one. 

From the time integrated data of ARGUS and CLEO, the time-dependent 
analysis of B-decays at the high energy colliders (LEP II, Tevatron, SLC) 
and the time-dependent analysis at the asymmetric 5-factories Belle and 
BaBar the following result was obtained [5]: 

x d = 0.755 ±0.015 . (11.6) 

By using the life time of B d -mesoias [5]: T Bd = [1.54(1) • 10~ 12 sec] -1 = 
[1.54(l)ps] _1 we get for the mass difference of B d mesons [5]: 

Am Bd = 0.489(8)ps -1 . (11.7) 

In Standard Model B d — B d transition occurs through the box diagram 
shown in Fig. 8. 

b u,c,t d 
< 



« u,c,t b 

Figure 8: Diagram responsible for B d — B d mixing 

Unlike the case of K° — K° transition the power of A is the same for u, c 
and t quarks inside a loop, so the diagram with t-quarks dominates. M 12 is 
the same as calculated in [21, 22] for the i^-meson case. 

From eq. (9.5) substituting B K f^V t 2 s by B Bd f B V t l and removing "Im" we 
obtain: 

M 12 = - G ^j l ' msm^BV^ 7(0 , (11.8) 
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where /(£) is the same function as that for J\~-mesons, eq. (9.5). i]b with the 
account of NLO corrections is [27]: 

r^ LO = 0.55 ±0.01 . (11.9) 

Ti2 is determined by the diagram analogous to that shown in Fig. 4 
(where s-quark should be substituted by 6-quark) but inside the loop c-quark 
can propagate as well (so, 4 diagrams altogether: uu, uc, cu, cc quarks in 
the inner lines). According to [28]: 

ri2 = G2 ^f>B [Vcb v: d (i + o{^)) + v ub v:f , (n.io) 

2 

where the term O(^) accounts for nonzero c-quark mass. 2 

b 

Using the unitarity of CKM matrix we get: 

ri 2 = ^y^ HW + oi^)v ch v: d f , (ii.ii) 

and the main term in Ti2 has the same phase as the main term in M12, eq. 
(11.8). That is why CPV in mixing of 5-mesons is suppressed by an extra 
factor (m c /mf,) 2 . Postponing the discussion of CPV in B — B mixing till the 
next section, from eq. (6.5)we obtain: 

M + - M_ - 1 -{Y + - r_) = 2[| M 12 I -- I T 12 |] , (11.12) 

and with the help of eq. (11.8) for the difference of masses of the two eigen- 
states we obtain: 

Am Bd = - ^yl ^ I V t lV td 2 I 7(0 , (11.13) 

and Am Bd is negative as well as in the kaon system. 

Comparing this theoretical formula with the experimental result (11.7) 
we obtain the bound on | Vu |, which forms the circle in Fig.l with the center 
at the point 77 = 0, p — 1: 

\Vu T= A 2 \ 6 [(l-pf + f] . (11.14) 

2 

2 The dependence of Fi2 on is more involved since the contributions of the diagrams 

b 

with uc and cc quarks depend differently on (m c /m,b) 2 . However this is not important for 
what follows. 
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The main uncertainties in this bound on rj, p plane are theoretical [5]: 

B B J 2 Bd = (1.3±0.1)(200±30MeV) 2 , (11.15) 

where the lattice results are used. 

Substituting the numbers in (11.13) using (11.7) we get: 



(l"P) 2 + f 



= (0.85±0.15) 2 . (11.16) 



If in the diagram shown in Fig. 8 we substitute s-quark instead of ci-quark 
we will get B s — B s transition. Making straightforward replacements in eq. 
(11.13) we obtain: 



B B d P Bd V td 



Am Bs = Am Bd -^f | -£ Am Bj (1.2 ± 0.1) 2 x 



x A2[(1 _ l p? + t] = ( 4 5 ± 7)Am Bd = (22 ± A)ps~ 1 , (11.17) 

where for the hadron parameters we use a lattice result from [5] , the numerical 
values of p and rj from the general fit (see Conclusions) and the experimental 
result for Am^. Since the lifetimes of Bd - and B s -mesons are almost equal, 
from (11.6) and (11.17) we get: 

x Bs =34 ±5 , (11.18) 

which means very fast oscillations. That is why W Bs equals 1/2 with very 
high accuracy and one cannot extract x Bs from the time integrated measure- 
ments. Performed up to now searches were not sensitive enough to measure 
x Bs ; only the lower bound was obtained [5]: 

Am Bs > U.lps' 1 at 95% CL . (11.19) 

New Tevatron run should be able to find the value of Am Bs if it is close 
to Standard Model prediction, eq. (11.17). 

However even the lower bound (11.19) appears to be powerful enough to 
bound possible values of p and rj quite effectively. From (11.17) and (11.19) 
we get: 

(l-p) 2 + f< 1.1 at95%CL , (11.20) 
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which produces Am# s circle in Fig. 1. One can see that this circle is al- 
ready inside the outer circle from (11.16) (recalculated to 2a deviation). 
This progress is possible because in the ratio (B Ba f B J/(B Bd f Bd ) theoretical 
uncertainty diminishes. 

For the difference of the width of B d + and B d - we obtain: 



Ar 



r 



12 



G%B 



B d 



which is very small: 



16tt 



0.3% 



td 



(H-21) 



11.22) 



as opposite to .fT-meson case, where K$ and lifetimes differ strongly. 

In the S s -meson system a larger time difference is expected; substituting 
V ts instead of V td in (11.21) we obtain: 



Ar, 



10% 



(11.23) 



12 CPV in B° - B u mixing 



»o 



dg L — too small 



For a long time CPV in i^-mesons was observed only in K° — K° mixing. 
That is why it seems reasonable to start studying CPV in _B-mesons from 
their mixing: 



Q 






P 




\ 




12 



M- 



12 , 




12 



M- 



12 , 



= 1 



-Im ( — — 

2 VM 12 



dm- 



v cb v; d 



mi 



1-O(10~ 4 ) 



(12.1) 



v tb v td 

where the expressions for M 12 (eq. (11.8)) and r 12 (eq. (11.11)) were used. 
We see that CPV in B d — B d mixing is very small because t-quark is very 
heavy and is even smaller in B s — B s mixing. The experimental observation of 
B d — B d mixing comes from the detection of the same sign leptons produced 
in the semileptonic decays of B d — B d pair from T(4S f ) decay. Due to CPV 
in the mixing the number of events will differ from that of / + / + and this 



difference is proportional to 



10" 



So, one needs more than 10 
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semileptonic decays of both B and B to detect this effect (statistical error 
~ \fN). Taking into account that semileptonic branching ratio of B to euX 
or fivX is 20%, we see that about f0 10 T(4S') — > B°B° decays are needed to 
observe CPV in mixing according to the Standard Model. This is completely 
hopeless now taking into account that 5-factories have collected about f 8 
BB pairs produced in T(4S') decays. 

Another type of experiment is possible as well. Let us suppose that in 
some reactions when B + B° pair is produced B + is detected, and if B~B° is 
produced B~ is detected. In this way a flavor state is tagged: it is known 
what particle (-B or B°) was produced. If neutral B decays semileptonically 
one can look for CPV charge asymmetry: 

4l = T7~~~~~ ~ Z B "^ l+ = 0(1(T 4 ) , (12.2) 

while the experimental result is (in [5] it is denoted by acp)'- 

af L = -0.002 ± 0.009(stat) ± 0.008(syst) - (12.3) 

- two orders of magnitude weaker bound on af L than the theoretical predic- 
tion. 



13 CPV in interference of mixing and decay 
(Im£ ji 0) 

As soon as it became clear that CPV in B — B mixing is small theoreticians 
started to look for another way to find CPV in B decays. The evident 
alternative is the direct CPV. It is very small in fT-mesons because: a) 
the third generation almost decouples in K decays; b) due to AT = 1/2 
rule. Since in £>-meson decays all three quark generations are involved and 
there are many different final states one can hope to have large direct CPV 
there [29] - [32]. An evident drawback of this strategy: a branching ratio of 
5-meson decays into any particular exclusive hadronic mode is very small 
(just because there are many modes available), so large number of 5-meson 
decays is needed. The specially constructed asymmetric e + e~-factories Belle 
and BaBar collected about thirty million BB pairs produced in T(AS) decays 
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each and discovered in 2001 CPV in B(B) decays [14] (at the end of 2002 
statistics is almost 3 times larger in each experiment). 

The time evolution of states produced at t = as B° or B° is described 
by eqs. (11.1). 

It is convenient to present these formulae in a little bit different form: 



M_i_+M_ rt 

B°(t) >= e- 1 ^— 1 -^ 



. .Amt, n .q . .Amt, 
[<x*(-2-)\B°>+i±sm(—, 



B°(t) >= e 



M_ 



. M i +M_ 



/-E*r -P ■ 

2 sm 



Amt, 



q 2 
M + > 0, and we put r + = T 



B° > + cos 



Amt, 



5 U >1 



5° >], 



(13.1) 

T neglecting their small 



where Am 
difference. 

Let us consider a decay in some final state /. Introducing the decay 
amplitudes according to the following definitions: 



Af = A(B° - /) , A f = A(B - /) , 

Af = A(B° - /) , A/ = A(fi - /) , 
for the decay probabilities as functions of time we obtain: 



(t) 



(t) 



-rt 



| 2 [cos 2 (— — ) + 



A/ | 2 [cos 2 ( 



Amt, 



) + 



sm 



, . Amt , 
-' )-/m 



gA 7 



sm 



9 . Amt , 

—)-Im 



'qA£ 

Ms. 
Mr) 



(13.2) 



sin(Amt)] , 



sin(Amt)] 
(13.3) 



The difference of these two probabilities signal different types of CPV: 
the difference in the first two terms in brackets appears due to direct CPV; 
that in the last term - due to CPV in interference of B° — B° mixing and 
decays. 

Let / be a CP eigenstate: / = rjff, where rjf = +(— ) for CP even (odd) 
/. (Two examples of such decays: B° — > J/^>Ks(l) an d B° — > 7r + 7r~ are 
described by the quark diagrams shown in Fig. 9. The analogous diagrams 
describe B° decays in the same final states.) The following equalities can be 
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easily obtained: 3 



A f = VfAf , A f = r) f Af 



(13.4) 




a) 



d d 



b) 



Figure 9: Quark diagrams describing B° — > J/^Ks(l) ( a )) and B° — > n + n 
(b)) decays 

In the absence of CPV the expressions in brackets in (13.3) are equal and 
formulas (13.3) describe the exponential particle decay without oscillations. 
Taking CPV into account and neglecting a small deviation of | p/q | from 
one, for CPV asymmetry of the decays into CP eigenstate we obtain: 



Pso^f - PB°->f I A I 2 -1 , . , 2Im\ 
a cp(t) = p — p = I \ |2 i-i cos(Amt) + sm(Amt) = 

= -C f cos(Amt) + 5/ sin(Amt) , (13.5) 



where 



(13-6) 

(do not confuse with the parameter of CKM matrix), Cf = |~|^| 2 , .5/ = 
ifpq^ • 4 The nonzero value of Cf corresponds to direct CPV; it occurs when 

3 The first one follows from the following identity: < / — rjff \ AV \ B° >= 0, where 
AV is the interaction Hamiltonian responsible for the decay; the second one follows from 
the analogous identity for B°. 

i i ri ^ 



4 Belle collaboration uses Sf and Af = —Cf 
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more than one amplitude contribute to the decay. For extraction of CPV 
parameters (angles of unitarity triangle) in this case the knowledge of strong 
rescattering phases is necessary. The nonvanishing Sf describes CPV in 
interference of mixing and decay. It is nonzero even when there is only one 
decay amplitude, and |A| = 1. Such decays are of special interest since the 
extraction of CPV parameters becomes independent on poorly known strong 
phases of the final particles scattering. 

The decays of T(4S) resonance produced in e + e~ annihilation are a pow- 
erful source of B°B° pairs. A semileptonic decay of one of the B's tags 
"beauty" of the partner at the moment of decay thus making it possible to 
study CPV. However the time-integrated asymmetry is zero for decays were 
Cf is zero. This happens since we do not know which of the two S-mesons 
decays earlier, and asymmetry is proportional to: 



The asymmetric _B-factories provide possibility to measure the time-dependence: 
T(4S) moves in a laboratory system, and since energy release in T(4S) — > 
BB decay is very small both B and B move with the same velocity as the 
original T(4S'). This makes the resolution of B decay vertices possible unlike 
the case of T(AS) decay at rest, when B and B being non-relativistic decays 
at almost the same point (for the detailed description of the experiments 
see the lecture of A. Bondar in these proceedings) [33]. The implementa- 
tion of the time-dependent analysis for the search of CPV in 5-mesons was 
suggested in papers [34, 35]. 



14 B° d (B° d ) -+ J/VK S (l), sin 2/3 - straight lines 



Currently this is the only decay channel where both 5-factories observe sta- 
tistically significant CPV. 

The tree diagram contributing to this decay is shown in Fig. 9a). The 
product of the corresponding CKM matrix elements is: V* b V cs ~ A\ 2 . Also 
the penguin diagram b — > sg with the subsequent g — > cc decay contributes 
to the decay amplitude. Its contribution is proportional to: 



oo 




—oo 



P ~ V us V* b f(m u ) + V cs V; b f(m c ) + V ts V t lf(m t ) 
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= V us V: b (f(m u ) - f(m t )) + V cs V; b (f(m c ) - f{m t )) , (14.1) 

where function / describes the contribution of quark loop and we subtracted 
zero from the expression on the first line just as we did when calculating the 
box diagram in Section 8. The last term on the second line has the same 
weak phase as the tree amplitude, while the first term has a CKM factor 
V us V* h ~ A 4 (p — ii])A. Since (one-loop) penguin amplitude should be in any 
case smaller than the tree one we get that with 1% accuracy there is only one 
weak amplitude governing B®(B®) — > J/^Ks(l) decays. This is the reason 
why this mode is called a "gold-plated mode" - the accuracy of the theoretical 
prediction of the CP-asymmetry is very high, and BriB^ — > J/^/K°) « 10 -3 
is large enough to detect CPV. Substituting in Eq.(13.5) |A| = 1 we obtain: 

a C p{t) = ImXsm(AmAt) , (14.2) 

where At is the time difference between the semileptonic decay of one of 
5-mesons produced in T(4S*) decay and that of the second one to J/^K S (l)- 
Using the following equation: 

A f = n f Aj , (14.3) 

where rjf is CP parity of the final state, we obtain: 




a=i-i [ bo ^ k ^ = m (144) 



A B ^j/^ Ks{L) \pj Bd A B o^ J/q , Ks(L) 

The amplitude in the nominator contains K° production. To project it 
on Ks(l) we should use: 

Wk {q)K 

getting {q)K in denominator. The amplitude in the denominator contains 
K° production, and using: 

K° = (14.6) 

(P)K 

we obtain the factor (p)x i n the nominator. Collecting all the factors together 
and substituting CKM matrix elements for Aj/Af ratio we get: 
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Since in P decays J/\I/ and Ks(l) are produced in P-wave, f]s(L) = — (+) 
(CP of J/^f is +, that of Kg is + as well, and (—1)' = —1 comes from 
P-wave; CP oi K L is -). 

Substituting the expressions for {q/p)B d and {pjq)K in (14.4) and taking 
into account t\s{l) we obtain: 

V td V tb V c b V cs VcdV cs 

which is invariant under the phase rotation of any quark field. From eq. 
(3.20) and Fig. 2 we have: 

sxg(V t lV td ) = 2n-P , (14.9) 

and we finally obtain: 



a C p(t) 



= (-) sin(2/3) sin(AmAf) , (14.10) 

J /* K S(L) 



which is a simple prediction of Standard Model. In this way the measure- 
ment of this asymmetry at P-factories provides the value of the angle (3 of 
the unitarity triangle. The results of Belle and BaBar are consistent; their 
average is: 

sin 2(3 = 0.73 ± 0.05(stat) ± 0.035(syst) . (14.11) 

This result is based on the analysis of approximately 80 • 10 6 pairs of PP 
produced in T(4S') decays per collaboration [36]. As a final state not only 
J/^K S (l) were selected, but also the other states with hidden charm 
(^'Ks,r] c K s ,XcKs)- The value of |A| was also determined from the absence 
of cos (Am At) term in asymmetry: 

|A| = 0.95 ± 0.035(stat) ± 0.025(syst) (14.12) 

in accordance with Standard Model prediction. From eq. (14.11) we obtain 
2 possible solutions with angle 2f3 in the first or second quadrant. Both 
are shown in Fig. 1 by straight lines and the first one coincides nicely with 
Standard Model expectations. 

Returning to formula (14.8) let us note that the decay amplitudes and 
K° — K° mixing do not contain complex phases, that is why the only source 
of it in CP-asymmetry in P° J/^/K decays is P° — P° mixing: 
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thus the phase comes from V td , that is why the final expression (14.10) con- 
tains the angle 2(3 - the phase of V tc i/V t * d . 




15 Bd — ► 7T + 7r , sin 2a, Penguin versus tree, 
I A/A |# 1 

The pair 7r + 7r~ produced in Bd decay has positive CP: CP(tt7t)i = o = +1. 
The tree level diagram contributing to this decay is shown in Fig. 9b. Let us 
suppose for a moment that it dominates in the decay probability analogously 
to B — > J/^/K case. Since CKM matrix element V u b has a nonzero phase 
the CP-asymmetry should be different from that in J/^K decays. Let us 
calculate it: 

\-( q \ ^ - VtdV tb VubV ud _ -i(2/3+2 7 ) n r-,N 

{p) B A-^/v^r d - e • (15 ' 1) 

Im\ = sin 2a , 

&cp{t) — sin(2a) sin(Amt) 

where the triangle relation a+fi+j = n is used. So, the study of t-dependent 
CP asymmetry in T(AS) — > BB — > 7r + 7r~ decay would measure angle 
ct if tree diagram dominates in i? — > 7T7t decay. The penguin diagrams pro- 
ducing the transition b — > dp — > d-u-u also contribute to the 7r + 7r~ decay 
mode: 

p ~ WW + ^KUK) + W/K) = 

= K^J/KJ - /K)] + V cd V; b [f(m c ) - f{m t )\ , (15.2) 

and while the first term should be added to the tree diagram (CKM phase 
is the same), the second one has a different phase and is of the order of 
A 3 , just like the tree diagram. Naively one would expect that the penguin 
contribution should be nevertheless much suppressed. Feynman diagram 
calculation leads to the following damping factor: 

P/T ^ h py«0.04 . (15.3) 
1 12tt \mj y ' 



42 



However the comparison of B — > txtx and B — > 7rfT branching ratios 
demonstrates that the naive estimate is not valid and the penguin contribu- 
tions are very much enhanced. The amplitudes of B — > nn decays are of the 
order of: 

A B ^ ~ A 3 T + A 3 P , (15.4) 

(here A is the CKM parameter, A ss 0.22), while that for P> — > 7T.K" decays 
are: 

A b ^k ~ A 4 T + A 2 P , (15.5) 

and if estimate (15.3) were correct, the decay probabilities with .fT-meson 
production should be suppressed as A 2 ~ 5 • 1CT 2 . Experimentally a number 
of branchings of these two types of decays were measured, or upper bounds 
were established (different decay modes have different charges of B, K and 
7r). It appears that the decays to ixK are even more probable than to nn [5]. 
This can happen only for P/T > A; it follows from the experimental data 
that: 

(P/T) cxp ~ . (15.6) 

Taking into account the penguin contribution instead of eq. (15.1) we get 
[37]: 



A 



e 2ia 



1 + 


p 

T 


e «(<5+7) 


1 + 


P 

T 


e »(5-7) 



;is.7) 



where the part of the penguin contribution proportional to V u dV* b is included 
into the tree amplitude and 5 = 5p — St is the difference of strong phases of 
the penguin and tree amplitudes, a and 7 are the angles of unitarity triangle. 

The experimental data on CP asymmetries in B — > 7r + 7r~ decay are 
currently controversial. According to Belle 



= -1.23 ±0.41 ±0.08 



C nn = -0.77 ±0.27 ±0.08 , (15.8) 

where the first errors are statistical, the second - systematical. Both numbers 
are different from zero at the level of 3 sigma. Since C n7T is different from 
zero the penguin contributions are not negligible. Extracting the penguin 
amplitude from the B — > nK decay branching ratios with the help of flavor 
SU(3) and the tree amplitude from the factorization hypothesis applied to 
B — > neu decay the authors of [39] conclude that \P/T\ ~ 0.3. Using this 
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result and the triangular relation 7 = tt — (3 — a one can extract the values 
of a and 5 from (15.8). The result is [38]: 

78° < a < 152° (15.9) 

for (3 = 23° and 0.15 < \P/T\ < 0.45. 

BaBar Collaboration did not observe CP violation in B d — > it + it~ decay 
[40]: 

= 0.02 ± 0.34 ± 0.05 

C nn = -0.30 ± 0.25 ± 0.04 . (15.10) 

We can hope that with more data available and corresponding diminishing 
of statistical error the Belle-BaBar controversy on CPV in Bd — > 7r + 7r~ decay 
will be resolved. 



16 Angle 7 

Angle (3 of the unitarity triangle is already measured with good accuracy; 
angle a will be determined from CP asymmetry in Bd — > 7r + 7r~ decays with 
better accuracy when more statistics will become available. Evidently the 
next task is to measure angle 7, or the phase of V u j,. In Bd decays angle 
(3 enters the game through Bd — Bd mixing. To avoid it in order to single 
out the angle 7 let us turn to B s decays: CKM matrix elements V ts and V tb 
which participate in B s — B s mixing are real. 5 Selecting the decays where 
b — > u transition dominates and looking for CP asymmetry we can measure 
angle 7. B s decays to CP eigenstates which occur through b — > u transition 
would provide the necessary information: B s — > Ksp, Ksii . However the 
penguin transition b — > dg — > duu is proportional to A 3 as well as the tree 
decay b — > uud. That is why the extraction of 7 from these decays will be as 
involved as that of a from B d — > decays. 

The analogs of "golden" Bd — > J j^K mode in the case of B s are B s — > 
J/^(j), J/^r] 1 decays. In the framework of Standard Model the CP asymme- 
tries in these decays should be zero: CKM phase is absent in B s — B s mixing 
and in b — * ccs decay amplitude (CP (J/^frf') = +, J/^<p i s a mixture of CP 
odd and even states, but this fact only dilutes asymmetry). That is why the 

5 Another way to avoid Bd — Bd mixing is to look for direct CPV in B* 1 decays. 
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search of this asymmetry is very interesting: a nonzero result discovers New 
Physics (most probably in B s — B s mixing). 

An interesting strategy for the measurement of angle 7 was suggested 
in paper [41]: the study of the time-dependent decay asymmetries of B s — > 
D+K-, B s -> D+K~ as well as B s -> D~K+, B s -> ^JX* decays. The 
diagrams contributing to the first pair of decays are shown in Fig. 10. The 
final states are not CP eigenstates (unlike the previously discussed B d (B d ) 
decays) . 




a) 



b) 



Figure 10: Dominant diagrams contributing to B s — > D+K (a) and B s — > 
D+K' (b) decays 

Tagging a parent meson we can study the time-dependent asymmetry of 
B s — > D+K' and B s — > D+K' decays and extract from it the following 
quantity: 

\ _M (161) 

A D+K--[-\ A V * V ' ^ ib -^ 

see eqs. (13.6), (14.4), where amplitude Ai corresponds to Fig. 10a, and 
amplitude A 2 corresponds to Fig. 10b (CKM matrix elements are written 
separately). Analogously from the study of B s — > D~K+, B s — > D~K+ 
decays we extract: 

_M A 2 v ub v; s 
Xd7K+ -\p) Bs aK^ • (16 - 2) 
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Multiplying these two lambdas we get: 

,1/* !/*!/, 

(16.3) 



x v x - (lY V jy± - P -«7 

x A D7^+ - I I T/* V 1/ V* 



the unknown hadronic amplitudes Aj cancel, (q/p)B s is real and phase 7 will 
be measured - the "only" problem is to collect enough statistics of tagged 
B S (B S ) — > D+K~ {D~ K + ) decays and to measure its time-dependent asym- 
metries. 

Another way to determine 7 is through B d decays into DK final states 
[42], see also [43]. 

Study of B ± and Bd decays into nK final states also allows to determine 
angle 7 [44]; the present experimental uncertainties do not allow to draw 
definite conclusions. 

Other strategies to measure (or constrain) angle 7 can be found in liter- 
ature; for the list of references look at a recent review [45]. 

17 CPV inB^ (j)K s , K + K K S , r)'K s : pen- 
guin domination 

These decays are dominated by penguin diagrams b — > sg — > sss. The 
diagram with an intermediate w-quark is proportional to A 4 , while those 
with intermediate c- and t-quarks are proportional to A 2 . In this way the 
main part of the decay amplitude is free of CKM phase, just like in case of 
Bd — > J/^K decays. A nonzero phase which leads to time-dependent CP 
asymmetry comes from Bd — Bd transition: 

a C p(t) = -?7/sin(2/3) sm(AmAt) , (17.1) 

analogously to B d -> J/^K decays, eq. (14.10) ([46], see also [47]). <j)K s 
and fj'Ks final states are CP-odd, while in case of the decay into three kaons 
final state is a mixture of CP-even and odd states. According to [48] CP-even 
final states dominate. Here are Belle results [48]: 

Mode sin 2(3 

(pK s -0.73 ±0.64 ±0.22 ( 

K+K-Ks 0.49 ±0.43 ±0.11^ ^ ' 

rj'K s ±0.71 ±0.37 ±0.05 , 
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where the third error for the K + K~ Ks mode arises from uncertainty in the 
fraction of the CP-odd component. The value of sin 2(3 obtained from CP 
asymmetry in B d — > (j)K s decay deviates from the Standard Model expecta- 
tion (eq. (14.1)) by more than 2a. According to BaBar [49]: 

sin 2(3(<pK s ) = -0.19 ± 0.52 ± 0.09 , (17.3) 

and the deviation is more moderate. A number of New Physics contribu- 
tions to the penguin diagram were immediately suggested which explain the 
deviation from Standard Model in (pK s mode observed by Belle Collabora- 
tion. More statistics is needed to clarify the present situation: perfect SM 
description of CP asymmetry in B d — > J/^/K decays which occurs at the 
tree level and possible New Physics contribution to loop-induced B d — > <f)Ks 
decay. Belle and BaBar differ not only in S<j,K s , but in C^k s as well: 

CVx s (Belle) = 0.56 ± 0.41 ± 0.16 

C <f)Ks (BaBar) = -0.80 ± 0.38 ± 0.12 (17.4) 



18 Conclusions: CKM fit and future prospects 

Four parameters of CKM matrix are fitted from the following experimental 
data: \V ud \ = 0.9734 ± 0.0008 [5], \V US \ = 0.2196 ± 0.00261 [5], \V cd \ = 
0.224 ±0.016 [5], \V CS \ = 0.996 ± 0.013 [5], |V c6 | = 0.041 ± 0.002 [5], \V ub \ = 
0.0036 ± 0.0007 [5], sin 2(3 = 0.73 ± 0.06, \e K \ = (2.282 ± 0.017) • 10" 3 , 
AM Bd = 0.489 ± 0.008 ps -1 . Here are the results of the fit: 

A = 0.223 ±0.002 

A = 0.82 ±0.04 

fj = 0.32 ±0.04 (18.1) 

p = 0.24 ±0.08 

X 2 /n.d.o.f. = 7.8/5 . 

The good quality of the fit is caused partly by avoiding controversial data, 
likeeq. (17.2). 
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For the angles of a unitarity triangle we obtain: 

P = 23° ± 3° 

a = 103° ±9° (18.2) 
7 = 54° ± 9° 

Future prospects are: 

1. Diminishing the experimental errors in f3 from B — > J/^/K and in a 
from 5 — > 7r + 7r~ CPV asymmetries - the better accuracy in CKM 
parameters; 

2. Measurement of asymmetries in B — » decay with better accuracy 
- check of SM; 

3. Measurement of AM Bs through B s — 5 S oscillation - check of SM; 

4. Measurement of the angle 7 - check of Standard Model. 

The rare K- meson decay K + — > -n + vv and K° — > 7r° z/z/ are also sensitive 
to the values of CKM parameters p and i] which make the measurements 
of the corresponding decay probabilities interesting: the deviations from SM 
predictions would signal New Physics. 

In this way the study of flavor physics in pre-LHC era could put an end 
to thirty years success of Standard Model. 

I am grateful to L.B. Okun for fruitful discussions, to E.A. Andriash and 
G.G. Ovanesyan for performing CKM fit and drawing Fig. 1, to Ya.I. Azimov, 
M. Gronau, T. Hambye, A. Soni and M.B. Voloshin for useful remarks. I am 
grateful to E.A. Ilyina for her help in preparing the manuscript. 

References 

[1] Cabibbo, N., Phys. Rev. Lett. 10, 531 (1963). 

[2] Kobayashi, M. and Maskawa, K., Prog. Theor. Phys. 49, 652 (1973). 



48 



[3] Harari, H. and Leurer, M., Phys. Lett. 181B, 123 (1986); Fritzsch, H. 
and Plankl, J., Phys. Rev. D35, 1732 (1987); Botella, F.J. and Chau, 
L.-L., Phys. Lett. 168B 97 (1986); Terazawa, H., Prog. Theor. Phys. 
63, 1779 (1980). 

[4] Leutwyler, H. and Roos, M., Z. Phys. C25, 91 (1984). 

[5] PDG: Hagiwara, K. et al, Phys. Rev. D66, 010001 (2002). 

[6] Morehouse, P. et al., hep-ex/0304019 (2003). 

[7] Narodetskii, I.M., this volume of SHEP(2004). 

[8] Abele, H. et al., Phys. Rev. Lett. 88, 211801 (2002). 

[9] Lee, T.D., Yang, C.N., Phys. Rev. 104, 254 (1956). 

[10] Ioffe, B.L., Okun, L.B., Rudik, A.P, ZhETF 32, 396 (1957). 

[11] Landau, L.D., ZhETF 32, 405 (1957); Landau, L., Nucl. Phys. 3, 127 
(1957). 

[12] Christenson, J.H., Cronin, J.W., Fitch, V.L. and Turlay, R., Phys. Rev. 
Lett. 13, 138 (1964). 

[13] NA 48 Collaboration, Fanti, V. et al., Phys. Lett. 465B, 335 (1999); 
K TeV Collaboration, Alavi-Harati, A. et al., Phys. Rev. Lett. 83, 22 
(1999). 

[14] BaBar Collaboration, Aubert, B. et al., Phys. Rev. Lett. 87, 091801 
(2001); Belle Collaboration, Abe, K. et al., Phys. Rev. Lett. 87, 091802 
(2001). 

[15] Sakharov, A.D., Pisma v ZhETF 5, 32 (1967); JETP Lett. 5, 24 (1967). 

[16] Nikolaev, N.N., Proceedings of the First ITEP Winter School 3, 30 
(1973) (in Russian); 

Frere, J.M., Surveys in High Energy Physics 9, 203 (1996); 
Nakada, T., Surveys in High Energy Physics 9, 217 (1996); 
Bigi, I., Surveys in High Energy Physics 12, 269 (1998); 
Lee-Franzini, J., Franzini, P., Surveys in High Energy Physics 13, 1 



49 



(1998) ; www.itep.ru/, XXV ITEP Winter School of Physics; 
Schubert, K.R., Surveys in High Energy Physics 16, 1 (2001); 
Nakada, T., Surveys in High Energy Physics 17, 3 (2002), www.itep.ru/, 
5th Moscow International School of Physics; 

Pakhlov, P., Surveys in High Energy Physics 17, 37 (2002), 
www.itep.ru/, 5th Moscow International School of Physics (2002); 
Bigi, I.I., Khoze, V.A., Uraltsev, N.G., Sanda, A.I., CP Violation, p. 
175, ed. C. Jarlskog, World Scientific, Singapore (1989); 
Shabalin, E.P, Physics-Uspekhi 44, 895 (2001); 
Quinn, H., hep-ph/0111177, Lectures given at Trieste (2001); 
Ligeti, Z., hep-ph/0302031, Lectures given at 2002 SLAC Summer In- 
stitute; 

BaBar Physics Book, http://www.slac.stanford.edu/pubs/slacreports/ 
slac-r-504.html; 

Proceedings of the First Workshop on the CKM Unitarity Triangle, hep- 
ph/0304132 (2003). 

[17] Pontecorvo, B., ZhETF 34, 247 (1957); Sov. Phys. JETP 7, 172 (1958). 

[18] Vainshtein, A.I., Zakharov, V.I., Novikov, V.A., Shifman, M.A., Yad. 
Fiz. 23, 1024 (1976); Sov. J. Nucl. Phys. 23, 540 (1976). 

[19] Herrlich, S. and Nierste, U., Nucl. Phys. B419, 292 (1994); Phys. Rev. 
D52, 6505 (1995). 

[20] Hambye, T., Kdhler, G.O. and Soldan, PH., Eur. Phys. J. CIO, 271 

(1999) . 

[21] Vysotsky, M.I., Yad. Fiz. 31, 1535 (1980); Sov. J. Nucl. Phys. 31, 797 
(1980). 

[22] Inami, T. and Lim, C.S., Progr. Theor. Phys. 65, 297 (1981). 
[23] Veltman, M., Nucl. Phys. B123, 89 (1977). 

[24] Jamin, M., hep-ph/9911390 (1999); Eeg, J.O., hep-ph/00 10042 (2000). 
[25] Albrecht, H. et al., Phys. Lett. 192B, 245 (1987). 
[26] Albajar, C. et al., Phys. Lett. 186B, 247 (1987). 



50 



[27] Buras, A.J. et al., Nucl. Phys. B347, 491 (1990); 
Urbam, J. et al, Nucl. Phys. B523, 40 (1998). 

[28] Hagelin, J.S., Nucl. Phys. B193, 123 (1981). See also Bigi et al., [16], 
page 202. 

[29] Anselm, A.A., Azimov, Ya.L, Phys. Lett. 85B, 72 (1979). 

[30] Bander, M., Silverman, D., Soni, A., Phys. Rev. Lett. 43, 242 (1979). 

[31] Carter, A., Sanda, A., Phys. Rev. Lett. 45, 952 (1980). 

[32] Bigi, I.I., Sanda, A.I., Nucl. Phys. B193, 85 (1981). 

[33] Bondar, A., these proceedings. 

[34] Dunietz, L, Rosner, J., Phys. Rev. D34, 1404 (1986). 

[35] Azimov, Ya.L, Uraltsev, N.G. and Khoze, V.A., Yad. Fiz. 45, 1412 
(1987); Sov. J. Nucl. Phys. 45, 878 (1987). 

[36] Belle Collaboration, Abe, K. et al, Phys. Rev. D66, 071102 (2002); 
BaBar Collaboration, Aubert B., et al., Phys. Rev. Lett. 89, 201802 
(2002). 

[37] Gronau, M. and Rosner, J.L., Phys. Rev. D65, 093012 (2002). 

[38] Belle Collaboration, Abe, K. et al, hep-ex/0301032 (2003). 

[39] Gronau, M., Rosner, J.L., Phys. Rev. D65, 013004 (2002). 

[40] BaBar Collaboration, Aubert, B. et al., Phys. Rev. Lett. 89, 281802 
(2002). 

[41] Aleksan, R., Dunietz, I., Kayser, B., Z. Phys. C54, 653 (1992). 

[42] Gronau, M. and Wyler, D., Phys. Lett. 265B, 172 (1991); Gronau, M. 
and London, D., Phys. Lett. 253B, 483 (1991). 

[43] Atwood, D., Dunietz, I., Soni, A., Phys. Rev. Lett. 78, 3257 (1997); 
Atwood, D., Dunietz, I., Soni, A., Phys. Rev. D63, 036005 (2001); 
Atwood, D., Soni, A., hep-ph/0206045; 0212071 (2002). 



51 



[44] Fleischer, R. and Mannel, T., Phys. Rev. D66, 054009 (2002). 

[45] Fleischer, R., hep-ph/0305267 (2003). 

[46] London, D. and Peccei, R., Phys. Lett. 223B, 257 (1989). 

[47] London, D., Soni, A., Phys. Lett. 407B, 61 (1997). 

[48] Belle Collaboration, Abe, K. et al, Phys. Rev. D67, 031102 (2003). 

[49] BaBar Collaboration, Aubert, B. et al., hep-ex/0207070. 



52 



